Rabi Spectroscopy and Excitation Inhomogeneity in a ID Optical Lattice Clock 
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We investigate the influence of atomic motion on precision Rabi spectroscopy of ultracold fermionic 
atoms confined in a deep, one dimensional (ID) optical lattice. We analyze the spectral components 
of longitudinal sideband spectra and present a model to extract information about the transverse 
motion and sample temperature from their structure. Rabi spectroscopy of the clock transition 
itself is also influenced by atomic motion in the weakly confined transverse directions of the optical 
lattice. By deriving Rabi flopping and Rabi lineshapes of the carrier transition, we obtain a model to 
quantify trap state dependent excitation inhomogeneities. The inhomogeneously excited ultracold 
fermions become distinguishable, which allows s-wave collisions. We derive a detailed model of this 
process and explain observed density shift data in terms of a dynamic mean field shift of the clock 
transition. 
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I. INTRODUCTION 

Evaluations of optical frequency standards based on 
thousands of ultracold neutral ®^Sr atoms confined in one 
dimensional (ID) optical lattices now surpass [l], Q the 
best ^•^^Cs primary standards [1, 0, 0- The high preci- 
sion achieved with the Sr optical lattice clock relies on 
high quality factor optical spectroscopy J6|, enabled by 
confining the atoms in a magic wavelength optical lat- 
tice [7] and probing along the strong confinement axis in 
the Lamb-Dicke and resolved sideband regimes [1, [§| . No 
studies yet indicate a fundamental limit at the 10~^^ level 
of accuracy d [13, Reaching new levels of accu- 

racy demands that ever more subtle spectroscopic effects 
are taken into account. Recently, a density related clock 
frequency shift was measured [10] , even though the clock 
is based on fermionic ^^Sr where collisions at tempera- 
tures of a few /iK are suppressed by the Pauli exclusion 
principle. 

In this paper, we present a systematic experimental 
and theoretical investigation of the Rabi spectroscopy 
process used in the *^Sr optical lattice clock. This inves- 
tigation results in a more detailed model of the density- 
dependent clock frequency shift. These results are used 
to describe our recent experimental work [lo| . which has 
led to several different theoretical models of fermion clock 
frequency shifts [III [13] ■ 

The experimental setup is described in Sec. |TT1 We 
present a perturbative model of the ID lattice in Sec lIIIi 
derive sideband spectra in Sec. lIVI and show how temper- 
ature can be obtained from their shape. Furthermore, we 
discuss time-dependent Rabi flopping and spectroscopic 
lineshapes as well as their dependence on sample temper- 
ature and probe-induced inhomogeneities in Sec. [V] In 
particular, we identify the effect of transverse motion on 
the excitation process as the main source of inhomogene- 
ity in the system, causing atoms to lose their indistin- 
guishability. Accurate modeling of the process in Sec. IVII 
allows relating this loss of indistinguishability to a dy- 
namic clock frequency shift proportional to the atomic 



density. 



II. EXPERIMENTAL SETUP 

The ultracold atomic sample of ^''Sr is produced by 
standard laser cooling techniques and trapped at the 
antinodes of a vertically oriented ID optical lattice. The 
^^Sr ^5'o — '^Pq clock transition is subsequently inter- 
rogated with laser light propagating coUinear with the 
lattice axis. The geometry is shown schematically in 
Fig- da)- After the spectroscopic probe has redistributed 
atomic population between the ground \g) (^Sq) and ex- 
cited Je) {'^Pq) clock states, the populations are measured 
via fluorescence detection on the ^Sq — ^Pi transition, 
heating the sample out of the trap. By repeating many 
such measurements, data is aggregated while the spec- 
troscopy laser is scanned across the clock transition to 
acquire a spectrum, the probe pulse time is varied to ob- 
tain population dynamics, or the laser frequency is sta- 
bilized to the clock transition for clock operation. 

Details of the setup are described elsewhere [H, IT6| . 
here we summarize the important experimental parame- 
ters for reference. The atoms are cooled in a two stage 
magneto-optical trap (MOT) on the transitions indicated 
by solid arrows in Fig. [TJb). The first stage uses the 
strong (30 MHz) '^Sq - '^Pi transition at 461 nm. The 
second stage MOT uses dual-frequency narrow line cool- 
ing [l3,[il] on the ^So{F = 9/2)-^Pi{F = 9/2, 11/2) in- 
tercombination lines (7 kHz) at 689 nm. The optical lat- 
tice beam is overlapped with the second stage MOT and 
atoms are directly cooled into the lattice. The ID optical 
lattice is formed by two counterpropagating laser beams 
near the Stark cancellation wavelength A = 813.43 nm, 
where the differential polarizability of ^Sq and "^Pq clock 
states is zero Q- In this experiment, a gaussian beam 
tilted at a small angle with respect to gravity is fo- 
cused to a waist wq ~ 30 fim (1/e^ radius of inten- 
sity) and retroreflected by a spherical mirror with match- 
ing curvature. The optical lattice forms in the overlap 
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between the two beams. The typical lattice depth is 
Uq — 130 X hi/y-ec — 22 fiK, where ly^^c — h/{2m.\^) is 
the lattice recoil frequency. 




FIG. 1: (color online), (a) Experimental setup. A one dimen- 
sional optical lattice is oriented vertically along the z-axis to 
suppress tunneling. The lattice wavelength is 813.43 nm. The 
clock probe beam is aligned coUinear with the optical lattice, 
and the probe beam is transmitted by the mirror used to retro- 
reflect the lattice, (b) The optical transitions used to cool are 
^So - ^Pi at 461 nm, and ^^o - ^Pi at 689 nm. The spec- 
troscopic transition between the ground \g) and excited |e) 
clock states is the doubly forbidden intercombination transi- 
tion ^5o-^Po at 698 nm. (c) High quality factor spectroscopy 
performed on lattice trapped atoms gives rise to three dom- 
inant spectral features; a central carrier (green arrow) where 
the motional state |n) is conserved, accompanied by a red and 
blue sideband (red and blue arrows) where the motional state 
decreases and increases by one, respectively. 

The sample temperature is controlled by additional 
cooling (heating) resulting in a sample of ~ 4000 atoms 
at temperatures variable from 1 to 5 /xK. The atoms 
are Doppler cooled (heated) along the transverse direc- 
tions, and sideband cooling (heating) is performed along 
the lattice axis with light detuned from the ^So{F ~ 
9/2) - ^Pi{F =11/2) transition. 

Using optical pumping, the atoms are then spin- 
polarized in a weak magnetic field. By choosing the cor- 
rect polarization for the optical pumping beam, atoms 
are polarized into one of the two maximally polarized 
states with magnetic quantum number mi;- = ±9/2. Dur- 
ing the polarization step, the sideband cooling (heating) 
light is simultaneously applied. 

For carrier spectroscopy relevant to optical clock oper- 
ation we probe the — "^-Po clock transition along the 
strong confinement axis of the lattice in the Lamb-Dicke 
regime and the resolved sideband limit [1, Q with a sub- 
Hertz cavity-stabilized diode laser ^1^. Here, the trap 
frequency of 70-80 kHz is much larger than the probe 
recoil frequency of 4.7 kHz and the natural transition 
linewidth of 1 mHz [2l|, This ensures a clock 

interrogation which is highly insensitive to Doppler and 
recoil effects. The duration of a clock spectroscopy pulse 
is 80 ms and the probe beam has a ^--^5 times larger waist 



than the lattice beam to minimize inhomogeneities in- 
troduced by the gaussian probe beam profile. To probe 
time-domain Rabi flopping, the pulse time is varied and 
the motional sidebands are probed by detuning the clock 
laser by the motional trap frequency. 

After applying the spectroscopy pulse, atoms remain- 
ing in the ground state are detected by measuring fluo- 
rescence on the ^Sq — ^Pi transition. This detection pulse 
heats all ^5*0 atoms out of the trap. The population in the 
•^Po state is then measured by first pumping the atoms 
back to the ^Sq state via ^Pi and by again measuring 
fluorescence on the 461 nm transition. Combining data 
from these two measurements results in a normalized ex- 
citation fraction. The optical lattice is reloaded with a 
new sample for each new measurement with a ~ 1 s cycle 
time. 



III. SPECTROSCOPY IN A 
ONE-DIMENSIONAL OPTICAL LATTICE 

Neglecting optical aberrations in the transverse beam 
profile, the resulting trapping potential near the waist 
can be written as 

U{z, r) ^ -Uo cos2(/fcz) e-2'-V'"o' + ^ngz, (1) 

where r = y'x^ + y^ designates the transverse distance 
from the lattice axis, z is the longitudinal coordinate, 
k = 2tt/X is the lattice wavenumber, g ~ 9.81 m/s^ is 
the gravitational acceleration, and m is the mass of ^''Sr. 
The resulting trap is a nearly vertical stack of flat ellip- 
soids ("pancakes") with an aspect ratio given by lattice 
wavelength and beam waist. 

The gravitational energy shift between neighboring 
pancakes (~ 1 kHz) breaks the translational symmetry 
of the potential and for our trap depths intersite tun- 
neling becomes strongly suppressed 23]. The resulting 
separation of a lattice band into isolated sites is called a 
Wannier-Stark ladder [^Jj. Each site has several states 
whose energies are determined by the gravity-free lattice 
structure, but their energy widths decrease according to 
the residual Landau-Zener tunneling rates to neighbor- 
ing sites. For our trap depths, temperatures of a few /iK, 
and low site populations, this complex problem simplifies 
drastically: we can consider single particles in isolated 
sites, neglect intersite tunneling, and use thermal aver- 
aging to evaluate the relevant spectroscopic parameters. 

However, even the single site potential is not separa- 
ble into independent coordinates and it has no analyti- 
cal solutions. We approximate the longitudinal potential 
in a site as a one-dimcnsional harmonic oscillator in z 
with a quartic distortion. The shallow transverse poten- 
tial given by the gaussian lattice beam profile is approxi- 
mated as a two-dimensional harmonic oscillator in r. We 
also include the first order coupling term between the 
longitudinal and transverse degrees of freedom and find 

[/(z, r) c Uoi-l + k^z^ + - ^z^ - ^zV). (2) 
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Treating the quartic distortion and the couphng term 
in first order perturbation theory for harmonic oscillator 
states \n) = \nxTny,nz) gives an energy spectrum 



En/h ~Vz{nz + 2^ + ^r{nx + ny + 1) 
- v^cJ^{nx +ny + l)(n2 + 



(3) 



In the following sections, we will model the most 
prominent spectroscopic features in a scan across the 
carrier and motional sideband resonances and compare 
these models to experimental data to extract informa- 
tion about sample temperatures. The temperature as 
well as an effective probe laser misalignment will then 
be used to quantify the degree of inhomogeneity in the 
system, leading to a model of coUisional frequency shifts 
for ultracold fermions in a ID optical lattice. 



where we identify the longitudinal and transverse trap 
frequencies from the harmonic approximation to the po- 
tential as 
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(4) 
(5) 



These relations allow determination of the trap depth 
and beam waist by measuring the trap frequencies. The 
number of states in the trap is approximately given by 



NzN^, with 



Un 



Nr 



N, — . 



(6) 
(7) 



Typical longitudinal and transverse trap frequencies in 
our experiment are 80 kHz and 450 Hz, respectively, so 
that Nz {Nr) ~ 6 (1000). Due to the quartic distor- 
tion by the sinusoidal lattice potential, the longitudinal 
energy gap 

^Vz~ Vrcc{nz + 1) - Vvcc — {nx + ny + 1), 

^z 

determines the measured longitudinal trap frequency as 
approximately Vz — J^rec instead of Uz- In contrast, the 
r^z^ coupling term in Eqn. [2] has a more subtle effect. 
Radially oscillating atoms explore trap regions with dif- 
ferent longitudinal trap frequencies and their response to 
the spectroscopy laser changes. 

Vibrational laser spectroscopy of trapped atoms is a 
well described topic The Rabi frequency associated 

with a traveling wave probe of wave vector kp between 
initial and final trap states \ni) and |ny) is 



f7o(n/|e' 



(9) 



where the free-space Rabi frequency ilo is given by 
the dipole matrix element between the electronic states. 
Thus atoms in different motional states will respond dif- 
ferently to the same probe field. The resulting spread in 
Rabi frequencies leads to inhomogeneities in the system 
for any non-zero temperature. 



IV. SIDEBAND SPECTROSCOPY 

The motional sideband spectrum of lattice-trapped 
atoms exhibits features similar to the well known spec- 
tra of harmonically trapped ions Q. Each data point 
in Fig. [5] was obtained by measuring the excited state 
fraction after choosing a probe laser detuning in the 
vicinity of the clock transition and applying an 80 ms 
Rabi pulse. The traces show a narrow carrier transition 
that is (mostly) free of motional effects, accompanied by 
motional sidebands at the longitudinal trap frequency 
[vz — 80 kHz). The red-detuned sideband (n^ — > — 1) 
is suppressed with respect to the blue-detuned sideband 
(jiz — > Uz + 1), indicative of the temperature along the 
strong confinement axis. In contrast to ion trap exper- 
iments, the sidebands are smeared out towards the car- 
rier [2^. This skewing is due to the coupling between 
the longitudinal and transverse degrees of freedom which 
makes the longitudinal transition frequency depend on 
the transverse motional state. The resulting line shape 
carries information about the potential and the longitu- 
dinal and transverse temperatures of the trapped atoms. 
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FIG. 2: (color online). Sideband spectra with longitudinal 
and transverse cooling (circles) and without in-trap cooling 
(open squares) by scanning the spectroscopy laser across the 
electronic transition with a Rabi pulse time of 80 ms and 
Rabi frequencies on the order of 1 kHz. The solid lines are 
combinations of fits to the carrier and sidebands. 

The clock transition's natural linewidth is on the or- 
der of mHz corresponding to a metastable lifetime on 
the order of 150 s [l^]. At large probe powers, satu- 
ration broadening can increase the carrier linewidth to 
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~1 kHz. At typical Rabi probe times of 80 ms, the pop- 
ulation dynamics have dephased to an equilibrium (see 
Sec.|V|, allowing a time-independent model of the side- 
band shape. In the following, we will focus on the first 
blue longitudinal sideband. 

The blue longitudinal sideband is produced by atoms 
undergoing the clock transition along with a motional 
transition \ni) = \nx,ny,n^) |n/) = \nx,ny,nz + I) . 
The probe laser detuning 5 exciting this transition is 
given by the energy difference in Eqn. [8] and depends 
on the transverse motional state. The relative blue side- 
band amplitude for a given probe detuning is thus given 
by the transverse motional distribution. By assuming a 
thermal distribution among the transverse trap states, we 
can relate the shape of the shallow sideband edge (fac- 
ing the carrier) to the transverse temperature, as shown 
in Appendix [X] We find an approximate line shape as 
a function of detuning 5 from the carrier (valid for the 
shallow sideband edge): 



blue 



{5) 



a 



X e[7(n, 



-q[1-5/7(»^)] 



(10) 
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with a 

linz) 



j^, uncoupled longitudinal energy gap 
(uz + 1), and the Heaviside function 



fe 

^z ^rcc 

G. The sideband shape in Eq. [TOlis approximately area- 
normalized and its main feature is an exponential slope 
falling towards the carrier, produced by the changing 
longitudinal energy gap with transverse excursion (the 
third term in Eqn. [8]). This exponential is multiplied 
by a linear term rising towards the carrier and vanish- 
ing at 6 = 7(^2), which arises from the degenerate two- 
dimensional transverse confinement. These two features 
dominate the shape of the sideband and capture the shal- 
low slope towards the carrier. Equation [TOl neglects the 
transition between the shallow edge and the underlying 
power-broadened Lorentzian that determines the sharp 
edge. 

The final, thermally averaged, blue sideband absorp- 
tion cross section is then given by (see Appendix E]) 



O-bluc(^) OC 



n-=0 



'bk 



(11) 



as a Boltzmann-weighted superposition of single (longi- 
tudinal) state sidebands shifted by the anharmonicity of 
the longitudinal trap. Here En, is the energy of longi- 
tudinal state Uz, neglecting the r^z^ coupling term. In 
addition, each — > + 1 sideband is smeared out to- 
wards the carrier by the coupling term between the lon- 
gitudinal and transverse traps. Each component's base 
Lorentzian would additionally be broadened by the life- 
time of the corresponding Wannier-Stark state. However, 
for the relevant states and trap depths, this broadening 
is much smaller than the power-broadened width and can 
be ignored. 



Regardless of the details of the component lineshapes 
in Eqn. 111! the only difference between the red and 
blue sidebands should be that the Boltzmann weights 
are shifted according to Uz f—^ + 1 since the particle 
starts in the higher-lying motional state. There is no con- 
tribution from the longitudinal ground state to the red 
sideband cross section and the ratio of integrated side- 
band absorption cross sections obeys 



^total 



-total 
^bluo 



(12) 



which can be solved numerically for Tz after determining 
Vz from the sharp sideband edge on the far side of the 
carrier, resulting in a lineshape-independent measure of 
the longitudinal temperature. 

Ideally, the transverse temperature Tr would be deter- 
mined with the same method by measuring the transverse 
sidebands at S — Vr, but their amplitude is suppressed 
since the probe is carefully aligned with the lattice axis. 
The transverse sidebands can be measured by misalign- 
ing the probe beam. An independent measurement of 
the trap frequency can be obtained by parametric heat- 
ing. Both methods confirm a typical transverse trap fre- 
quency of Vr — 450 Hz. A combination of a longitudinal 
sideband spectrum with a transverse trap frequency mea- 
surement allows a complete characterization of the trap 
parameters including depth Uq and waist wq via Eqn. |4l 
This calibration, together with a fit of the longitudinal 
sideband shape with Eqn. 1111 then lets us extract infor- 
mation about Tr from the same data that gives Vz and 
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FIG. 3: (color online), (a) Blue sideband as a function of 
probe time. The solid curves are fits of Eqn. Illl to the data 
for fixed Tz — 3 /xK determined from the cross section ratio 
in Eqn. 1121 (b) Apparent transverse temperature as deter- 
mined from the fits in (a) as a function of probe time on a 
double- logarithmic plot. The apparent transverse tempera- 
ture increases with probe time, although unity-aspect TOF 
data indicates Tr = Tz (dashed line). The solid line con- 
nects the first and last point, giving an increase in apparent 
transverse temperature as ~1.5 fiK per decade of probe time. 

Spectroscopic lineshapes of the blue sideband are dis- 
played in Fig.[31[a), along with fits of Eqn. [11] for varying 
Rabi pulse time in the range of 0.5-32 ms. The longitu- 
dinal temperature was determined by the ratio of red 
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and blue sideband areas (Eqn. [T2|) as Tj = 3 /iK and 
time-of-flight (TOF) pictures show a unity aspect ratio, 
indicating a well thermalized sample with — Tr- Fig- 
ure[2fb) shows the transverse temperature resulting from 
the fits in (a) and the dashed line indicates T^; the solid 
line connects the first and last points. 

The data show an apparent increase in transverse tem- 
perature with increasing probe time. A similar mea- 
surement with constant pulse time and varying probe 
power shows no such dependence. In combination with 
the TOF measurements, we conclude that the transverse 
temperature is not influenced by the probe laser, and that 
the magnitude of the Rabi frequency does not influence 
the sideband shape except by broadening the underly- 
ing Lorentzian. The time-dependent data suggest that 
the transverse atomic motion lets more atoms interact 
with the probe laser than suggested by the steady-state 
eigenenergy picture that gives rise to Eqn. [TOl With 
probe pulse times on the order of the inverse trap fre- 
quency, a snapshot of the atomic distribution can be ob- 
tained and the sideband shapes derived above become 
more applicable. 



V. EFFECT OF RADIAL MOTION ON THE 
CARRIER TRANSITION 

The longitudinal sideband shape is strongly modified 
by the anharmonic corrections to the potential. In con- 
trast we expect the carrier to be insensitive to motional 
effects since carrier spectroscopy can be done at line 
widths of ^2 Hz, corresponding to spectroscopic qual- 
ity factors of 2 X 10^"' 6], which is well in the resolved 
sideband regime. 

The carrier Rabi frequency for a harmonically trapped 
particle in state \n), probed with a spectroscopy beam 
of wave vector kp — 27r/Apfcp for probe wavelength Xp is 
given by 

nr. = no{n\e'''--^\n)=no [] e^"- (r,^), (13) 

i—x.y,z 

determined by the Lamb-Dicke parame ters rii = 
kpai/\^, with oscillator length 27rai — y/h/ [mvi) for 
trap axis i and Laguerre polynomial ■ In this way, 
the Rabi frequency is determined by the spatial extent 
of the motional wave function with respect to the probe 
wavelength. For a probe that is perfectly aligned with the 
lattice axis, the only non-zero Lamb-Dicke parameter will 
be 77z- In this case, -ql = v^^Jv^, where v^^^ ~ 4.7 kHz is 
the probe recoil frequency. 

We ignore corrections to Eqn. [13] from perturbed har- 
monic oscillator states, since both the and r^z^ terms 
in Eqn. [2] produce wave function coefficients suppressed 
by t'roc/(2i'z). Although the r^z^ coupling modifies r\z^ 
the correction to its value scales as x ^v-c^^v^jv;, and is 
negligible for the relevant temperatures. In the following, 
we will assume harmonic oscillator states and energies 



and derive carrier lineshapes and Rabi flopping accord- 
ingly 

Although we can ignore the trap shape except for the 
harmonic confinement, the Rabi frequency expression 
still assumes that the particle is probed by a plane wave 
without transverse profile. Each wave vector component 
of the probe beam contributes to the Rabi frequency ac- 
cording to its Fourier coefhcient, allowing us to estimate 
the effect of a shaped probe beam. Over the extent of 
a trap site the probe beam shape can be approximated 
as a plane wave along the mean probe direction with a 
transverse shape function describing its intensity profile. 
Although the probe beam is carefully aligned with the 
lattice axis, residual misalignment or aberrations might 
produce a net mismatch between lattice axis k = kz 
and probe axis. Since the transverse trap is isotropic, 
we choose a small net misalignment angle Ad along x, 
such that kp ~ kp{z + A6x). The transverse extent of 
the probe beam is large and can be approximated by a 
cylindrically symmetric function (with respect to the net 
probe direction) of waist Wp 3> wq. The corresponding 
Rabi frequency is 

- r!o(n|e"=''- [1 + 0{al/wl)] \n) 
-noinle'^^-^ln) (14) 

Thus the transverse shape of a cylindrically symmetric 
probe beam with large cross section cannot influence the 
spectroscopy {a^/Wp ~ 10~^ for our experiment). The 
Lamb-Dicke parameters are 



77, = l/\p^h/{2mM,) (15) 
= Ae/Xpyyh/{2mj^r) (16) 

Since i>r ^ ^z, even a small amount of effective misalign- 
ment will cause significant transverse contributions to the 
carrier lineshape. For our trap frequencies, 77, ~ 0.24, 
and rjx ~ AO x 3.2. In conclusion, the main correction 
to the carrier Rabi frequency comes from an effective 
misalignment angle in the transverse direction, while the 
broad transverse profile can be ignored. This misalign- 
ment introduces information about the transverse mo- 
tional state distribution and thus the transverse temper- 
ature into the carrier line shape. We will present experi- 
mental evidence for a typical effective misalignment angle 
Ae ~ 10 mrad. 

The carrier line shape for Rabi spectroscopy can be 
understood as follows. We can neglect spontaneous emis- 
sion and probe laser decoherence even at Rabi spec- 
troscopy times approaching 1 s Q. The lattice lifetime 
has also been measured as ~ 1 s and does not introduce 
decoherence on the typical Rabi pulse times of 80 ms. In 
this regime, we can neglect any decoherence rates in the 
system and describe the population dynamics in a fully 
coherent way. The excited state probability for motional 
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state |n), detuning S and pulse time t is 



Pe{n,d,t) 



"-"^ sin^ 



52 



(17) 



The ensemble-averaged excited state population is then 
given by 

Pe(<5, - E ^^r)qn^ {Tz)Pe{n, S, t) , (18) 



for normalized Boltzmann- weights qn^ [q-n^ ) correspond- 
ing to transverse (longitudinal) temperature (T^): 



^n, = (1 - Zx)zl 

qn, = (1 - 2^)2:^ 



= exp[-/ii/r/(fc_BTV)] (19) 
EE eM-hvz/{kBT,)] (20) 



We are interested in two scenarios: (a) Change the pulse 
time t and measure Rabi flopping at zero detuning, (b) 
Change detuning 5 by scanning the spectroscopy laser 
across the carrier resonance and measure the Rabi line 
shape. Both cases are covered by Eqn.[T8]and we expect a 
response as a coherent superposition of slightly different 
Rabi frequencies. 

For low temperatures, we only have to consider atoms 
in the lowest longitudinal state 71,2 = and for small 
misalignments, the Rabi frequency can be expanded in 
?7^. These two simplifications allow finding an analytical 
expression for case (a). Here, the excited state population 
can be approximated as 



00 

_ 1 l_ Zx Zx cos[0(l --ql)]- COS (f) 
~ 2 2 



(21) 



1 + zl~ 2zx cos{<j)r]l) 



with = 27rtfioe"'''/2e-'?^/2. Equation ED reduces to 
Pe{t) = sin^(/)/2 for zero misalignment and exhibits de- 
phasing to 1/2 by introducing different frequency com- 
ponents via r/2(A0) and amplifying their contribution by 
increasing Zx{Tr). For case (b), the sum can be evaluated 
in a similar manner. The resulting expression does not 
provide much further insight and we omit it here. 

Experimentally, the excited state fraction is obtained 
by measuring the number of atoms left in the ground 
state after the spectroscopy pulse, then repumping the 
excited atoms back to the ground state and measuring 
their number. The repumping back to the ground state 
has efficiency /3 < 1 and can vary from day to day. To 
obtain an excited state probability independent of a fluc- 
tuating overall number of atoms, the fraction of excited 
state over excited plus ground state counts is used. This 
normalization results in a Rabi flopping trace that de- 
phases to /?/(/? -I- 1) instead of 1/2. Note that on the 
experimental timescale, the dephasing is a coherent pro- 
cess dominated by an inhomogeneous distribution of Rabi 
frequencies via effective probe misalignment and temper- 
ature. 



We investigated the effect of inhomogeneous excita- 
tion by mapping out the Rabi flopping and lineshapes 
under different experimental conditions such as sample 
temperature, misalignment angle and probe laser inten- 
sity. In Fig. m^a) we show data for the excited state 
fraction as a function of probe pulse time for two differ- 
ent temperatures. For the hot sample, Rabi oscillations 
quickly decay after only 2-3 cycles. For the cold sample, 
the Rabi oscillations have a significantly higher visibility 
and are observable for about 10 cycles. This behavior 
was reproduced for a wide range of probe laser intensi- 
ties. In Fig. Sfb), the Rabi oscillations dephase faster 
when the effective misalignment angle is increased 
by misaligning the probe with respect to the lattice axis. 
From fits of Eqn.[2T]we determine a misalignment angle 
of A6' — 10 nirad for the well-aligned case relevant to 
clock operation. 

Each data point in Fig. 2] is determined by setting a 
specific probe pulse time and then scanning the probe 
laser detuning across the carrier transition to obtain a 
spectrum such as the one shown in Fig. [5] for a probe 
time of 1.7 ms. The maximum excited state fraction at 
zero laser detuning (~ 0.65 here) would then give the cor- 
responding data point in Fig. |3] for the respective probe 
time. 

For small A0, the data shown in Figs. [4] and [5] are 
well reproduced by the simplified model in Eqn. 1211 This 
model allows us to determine the transverse temperature 
as 1 /J,K for the cooled case, and 3 /iK for the uncooled 
case, with a precision of about 0.2 /iK, agreeing well with 
the corresponding longitudinal temperatures. In the mis- 
alignment model, both angle and transverse temperature 
produce a similar effect so that their covariance is sig- 
nificant. The transverse temperature data has been con- 
firmed by the unity aspect ratio of the atomic cloud in 
time-of-flight expansion. Accurate determination of the 
effective misalignment is difficult, introducing larger un- 
certainty in the radial temperature measurement. How- 
ever, fitting both Rabi flopping at zero detuning and Rabi 
lineshapes across a wide range of parameters gives con- 
sistent results. 

We characterize the amount of inhomogeneity in the 
system by the ratio of RMS spread in Rabi frequency 
Ari over the site's mean Rabi frequency il, given by 

= E q{nx)q{nz)^n^,n, (22) 

Tlx -iTiz 

An^ = E q{nx)q{n,)nl^^,,^ - n\ (23) 

In Fig. [HI we show how the ratio Afi/fi changes with 
sample temperature and misalignment angle. For the 
temperatures and effective misalignment angles, we find 
typical values of A$7/f2 in the range of 0.05-0.4. 

An additional source of Rabi frequency inhomogeneity 
between different lattice sites is the imperfect transmis- 
sivity of the lattice reflector for the probe wavelength, in- 
troducing a standing wave component in the probe beam. 
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FIG. 4: (color online). Effect of inhomogeneity on Rabi oscil- 
lations, (a) The excited state fraction as a function of pulse 
time is shown for sample temperatures of 1 /iK (circles) and 
3 fiK (open squares) with fits of Eqn. 1211 giving bare Rabi 
frequencies Qq = 59 Hz and 76 Hz, respectively. Both fits re- 
sult in a misalignment angle AO ~ 10 mrad. (b) The excited 
state fraction as a function of pulse time at 3 /iK is shown for 
increased misalignment angles and probe powers. The open 
circles show oscillations when the probe beam is slightly mis- 
aligned. The fit gives Qo = 55 Hz and AO ~ 10 mrad. The 
open squares show data for increased misalignment and the 
same probe power, giving Slo = 54 Hz and AO ~ 17 mrad. 
Note that the small misalignment approximation used to de- 
rive Eqn. [5T] starts to break down, resulting in a worse fit. 
Finally, the triangles show data for a large misalignment an- 
gle, taken with increased power. Although the fit looks much 
worse, the resulting AO ~ 40 mrad agrees well with a ge- 
ometrical estimate based on the experimental procedure for 
misaligning the beam. 
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FIG. 5: (color online). Excitation fraction as a function of 
probe laser detuning for a probe time of 1.7 ms. The solid 
curve is a fit of Eqn. [TS] with the temperature as the only free 
parameter, fixing at 10 mrad. The fit gives Tz = Tr = 
(2.1 ± 0.2) consistent with the sideband method and time- 
of-fiight expansion. 
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FIG. 6: (color online). Contour plot of the relative spread 
in Rabi frequency AQ/Q, as a function of misalignment an- 
gle AO and temperature Tz — Tr = T. The dashed contour 
lines indicate points where the mean Rabi frequency becomes 
smaller than the RMS spread, made possible by the negative 
values of the Laguerre polynomials in Eqn. [14] for large ra- 
dial quantum numbers. The experimentally relevant range is 
An/Cl = 0.05 - 0.4. 



The wavelength mismatch between the lattice and probe 
beams leads to a modulation of the residual probe stand- 
ing wave intensity between different lattice sites, and 
therefore different Rabi frequencies. However, this ef- 
fect is different from the excitation induced inhomogene- 
ity discussed above since it is homogeneous within each 
site. The mirror reflectivity at the probe wavelength was 
measured as 0.5%, leading to a Afl/Cl contribution of 
0.07. In the following section, we will use the inhomo- 
geneity information obtained from spectroscopy to es- 
timate inhomogeneity-induced clock frequency shifts re- 
lated to in-site Rabi frequency variation. The Rabi flop- 
ping data thus overestimate the in-site Rabi frequency 
variation slightly. 



VI. INHOMOGENEITY-INDUCED DENSITY 
SHIFTS 

For ^''Sr at /iK temperatures, interatomic collisions are 
suppressed due to Fermi statistics, which exclude colli- 
sions between identical fermions resulting from even par- 
tial waves. However, using ^''Sr a non-zero density shift 
to the transition frequency was recently measured [10| . 
At ultracold temperatures p-wave and higher order odd 
partial wave collisions are frozen out, and any remaining 
s-wave interactions should be suppressed. As discussed 
in the preceding sections, the trapping potential will re- 
sult in an inhomogeneous excitation process unless all 
atoms are in the same motional quantum state. Since 
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they are fermions that have been prepared in the same 
electronic and nuclear states, this is in principle impos- 
sible and there will always be an inhomogeneous excita- 
tion unless a lattice with at most one particle per site 
is realized. The inhomogeneous excitation changes the 
quantum statistics of the trapped atomic sample. This 
change in statistics leads to distinguishable fermions even 
for /xK temperatures and gives rise to s-wave collisions, 
and therefore a density shift. However, we will show that 
these small density-related frequency shifts can be mini- 
mized by choosing an average excitation fraction close to 
50%. 

The mean field energy shifts in the two state system 
consisting of ground and excited clock state \g) and |e) 
can be expressed as [l^ 

Ai?, - ^(Gg'aegpe + Gfjag.pg) (24) 
AE, - l!^(Gi2)a,,p, + Gg^aegPg), (25) 

(2) 

where G\j (cnj) is the two-body correlation function 
(scattering length) between states i and j, and pj is the 
partial density of state j in a trap site. The above expres- 
sion is applicable when the de Broglie wavelength of the 
scattering particles is large compared to their scattering 
length. 

The clock frequency shift can be calculated from the 
differential energy shift AE^ — AEg. For identical 

(2) 

fermions vanishes. Its bosonic equivalent equals two 
and for degenerate bosons it becomes unity. The clock 
frequency shift Av in the two-state fermionic system is 
given by 

hA,.= ^a,G(^J[pg-p,]. (26) 

We can estimate the magnitude of this shift by consid- 
ering the evolution of two representative atoms during 
the Rabi pulse. These two atoms will follow slightly dif- 
ferent trajectories on the Bloch sphere since their Rabi 
frequencies differ through the inhomogeneities derived in 
the previous sections, causing the energy difference be- 
tween the clock states to change during the clock pulse. 

Under the clock pulse, the electronic state of particle 
j evolves with the time-dependent Hamiltonian 

H,{t) = -'^{n^a, + [S- Ai^{t)]a,}, (27) 

where ax (<Tz) is the first (third) Pauli matrix, ilj is the 
Rabi frequency for particle j, and 6 is the probe laser 
detuning from the bare atomic resonance. We assume 
that all particles are prepared in the excited electronic 
state |e), and use a discretized version of the Schrodinger 
equation to propagate the wave function coefficients over 
a small time step At: 

\^j{t + At)) = {h - iAtHj{t)/h}\ijj{t)), (28) 



where I2 is the 2x2 identity matrix. After each time step 
At, the density shift Ai^{t) is recalculated for the updated 
wave functions |f/'j(i)) = o^jit)]^) + Pj{t)\g) (j ~ li2), 
according to Eqn. 1261 The antisymmetrizcd two-body 
correlation function between the two-level fermions 1 and 
2 is [3 

Gg) = 1 - \aiit)a;it) + Mmit)\^- (29) 

A more detailed many-body model relating G^^-' to tran- 
sitions to motional singlet and triplet states in the wave 
fimction is presented in reference 

If we consider a "mean particle" (particle 1) evolving 
with Rabi frequency Cl and a typical perturbing parti- 
cle (particle 2) evolving with Cl + Aft, we can define an 
approximate associated density shift as 

Ai^it) ^ Ai^o[l - \aia; + - 2\ai\% (30) 

where Avq = 2hpoaeg/m, with average atom density pQ. 
Here, the two-body correlation function is multiplied by 
the inversion of the mean particle following a trajectory 
on the Bloch sphere. Due to the background interac- 
tion with the other particles, a mean field energy builds 
up during the pulse, causing the detuning and thus the 
effective Rabi frequency to change dynamically. Many 
collisions occur during the Rabi time since the rate of 
collision attempts is given by the inverse of the trans- 
verse trap frequency l/i'r — 2 ms and the typical Rabi 
time for clock operation is 80 ms. This separation of time 
scales makes the mean-field treatment applicable. For 
pulse times <^ ^/i^r, the above description of the excita- 
tion process would have to be modified to include wave 
packets of vibrational states. For these short pulse times, 
no collisions can occur which can be explained by a local 
picture. The atoms are effectively confined to a volume 
determined by their velocity and the pulse length. If this 
volume is significantly smaller than the trap size, the lo- 
cal excitation process is very homogeneous. No collisions 
can occur since the atoms cannot travel far enough to 
encounter collision partners that have been excited in a 
slightly different way. 

The clock transition frequency is measured by locking 
the spectroscopy laser to points of equal height on the 
transition lineshape. To model the experimental proce- 
dure, we include the mean-field density shift as a time- 
dependent detuning in the wave function coefficients via 
Eqn. [28l We calculate the line shape for a pulse time 
corresponding to a 7r-pulse (on resonance for the mean 
particle) as a function of the inhomogeneity parameter 
Afl/ri, as shown in Fig. [7]Ja). The measured clock shift 
can be visualized by comparing two contour lines of equal 
excited state fraction pe- The spectroscopy laser is locked 
to their average detuning which changes with both the 
inhomogeneity and the pair of contours chosen. The fi- 
nal frequency offset is negative for pg < 0.52, vanishes 
at Pg ~ 0.52 and becomes positive for pg > 0.52. The 
resulting frequency shift (or locking error) with respect 
to the bare atomic transition in units of Avq is shown in 
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Fig. [T^b) for the experimentally relevant range of AJl/fj. 
The frequency shift is shown with a negative sign to recre- 
ate the behavior for a negative scattering length Oeg . 

In using Eqn. [301 and detecting only the mean particle 
excitation fraction, we have chosen a maximally asym- 
metric model where the representative particle is only 
perturbed by the background particles but does not act 
back on the perturbing particles. A maximally symmet- 
ric model can be derived by choosing representative Rabi 
frequencies f2 ± Af2, and using the mean of both excited 
state fractions in Eqn. [30] as well as in the detection of the 
line shape. In this two-particle model, there is no pre- 
ferred particle and the resulting lineshape exhibits loss 
of contrast as the inhomogeneity grows. We find that 
the density shift calculated from such lineshapes shows 
a very similar dependence on the ground state fraction 
and its magnitude is the same. 

To extract information about aeg from experimen- 
tal data for small inhomogeneity, we fit the numerical 
clock shift data with a with a simple polynomial. The 
frequency shift increases quadratically with Afi/fJ, is 
quadratic in pg and changes sign close to Pg = 0.5. Fit- 
ting the functional form 



A{pg - B){pg - c) {An/ny 



(31) 



to numerical data obtained after a 7r-pulse excitation, as 
in Fig.EIb) for AU/n e [0, 0.4], results in A = -0.410(3), 
B = 0.528(1), and C = -0.443(7) with RMS error 3 x 
lO^''. To a good approximation, we can thus use the 
formula 

Au^ ~ -Avq [AVL/nf QAl{pg - 0.53)(pg + 0.44) (32) 

to calculate the density shift for small inhomogeneity, 
measured by locking the spectroscopy laser to points of 
equal ground state fraction after a 7r-pulse (defined on 
resonance for the mean particle) for particles initially pre- 
pared in the excited electronic state. However, if the tt- 
pulse condition is allowed to vary by only 5%, the density 
shift curve changes as indicated by the shaded areas in 
FigEKb). Since we cannot exclude such variation over the 
course of a day, the zero crossing is not well defined and 
obtains an uncertainty according to how well the probe 
intensity is stabilized on long timescales. In the following, 
we will use the standard error given by fitting Eqn. [321 
to experimental data using weighted least-squares to in- 
dicate experimental uncertainty and add an estimated 
relative amplitude error of 30% in quadrature to account 
for the TT-pulse condition uncertainty. 

The density shift data from Campbell et al. [10] is 
shown in Fig. [8] along with fits based on Eqn.[32l adjust- 
ing the amplitude of the polynomial in pg as the only fit 
parameter. For temperatures Ti = 1 /iK and T2 — i /iK, 
we find fit amplitudes of Av{Ti) = -2.3 ±0.7 Hz and 
Au{T2) = -12.1 ±0.9 Hz, respectively. Although we 
can exclude p-wave scattering contributions (To| , the cor- 
responding values of aeg are larger than the unitarity 
limit given by a^g = Xt/'^t^, with thermal wavelength 
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FIG. 7: (color online), (a) Contour plot of the ground state 
fraction pg = 1 — |aip for the mean particle as a function of 
detuning in units of f2 and inhomogeneity parameter AQ./Q, 
after a vr-pulse. The solid curves are contours of equal pg, 
visualizing the locking points of the spectroscopy laser. Note 
that the contours for locking close to the Full Width Half 
Maximum (FWHM) are insensitive to the inhomogeneity pa- 
rameter, (b) The solid lines show the clock frequency shift 
(or locking error) with respect to the bare atomic transition 
in units of Aj/q = 2haegpo/m as a function of final ground 
state fraction pg = 1 — |cn|^ after the pulse for the relevant 
range of AQ,/Q.. A negative sign was included in the shift 
to recreate the experimental behavior for a negative scatter- 
ing length. The shift crosses zero close to pg =0.5 and varies 
quadratically both with inhomogeneity parameter and ground 
state fraction. The shaded areas indicate the change in the 
solid curve when the vr-pulse condition is allowed to vary by 
5%. 



At = h/ ■\/27rmfcB(r + T^p) and zero-point temperature 
Tzp ~ 3.5 /iK corresponding to the ground state energy 
of the trapping potential. On the other hand, in the uni- 
tarity limit, the ratio of both shift coefficients scales as 



(A17/f])i 



{An/il)2 



Ti + Tzp 



0.14, 



(33) 



which agrees with the ratio of measured shift coeffi- 
cients _Ai/(ri)/Az/(T2) = 0.2(1) for {An/n)i = 0.05 and 
(Arj/f2)2 = 0.14 at Ae = 10 mrad. 

The current two-particle model describes the qualita- 
tive behavior of the experimental data well and the scal- 
ing factors derived agree with the unitarity limit, but 
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FIG. 8: (color online). Clock frequency shift for 1 jjK (circles) 
and 3 fiK (open squares), normalized to the typical operating 
density po = 1 x 10^^ cm~^. The solid lines are fits of Ean.l32l 
to the data with the amplitude as the only fit parameter. 



cannot explain the large density shift quantitatively. This 
discrepancy probably stems from applying Eqn. [26] in a 
dynamic context. Even though there are many collisions 
during a typical Rabi pulse time, Eqn. [26] describes the 
mean field differential energy shift in a two-level system 
in steady state and cannot model the lineshape evolu- 
tion during the Rabi pulse accurately. A full many- 
body model is expected to improve our understanding 
of the density shift considerably JJJ . Even with a better 
model, extracting reliable information about the scatter- 
ing lengths from density shift data remains a challenging 
problem. 



VII. CONCLUSION 

In this paper we have investigated the inhomogeneous 
excitation introduced by the transverse degrees of free- 
dom for Rabi spectroscopy in a ID optical lattice. An 
analytical model for how the longitudinal sidebands are 
influenced by the finite temperature of the atomic sample 
is developed, and is used to allow us to extract temper- 
atures from sideband scans. Furthermore, by modeling 
and fitting Rabi oscillations and lineshapes the degree 
of inhomogeneity introduced by transverse temperature 
and an effective probe beam misalignment can be mea- 
sured. The inhomogeneity directly affects the spectro- 
scopic process and causes otherwise identical ultracold 
fermions to become distinguishable during the excitation 
process. This effect is modelled by a time dependent 
two-particle correlation function giving rise to density- 
dependent clock frequency shifts. By determining the 
temperature and inhomogeneity of the atomic sample as 
described here, and by controlling the excitation fraction, 
this density shift can now be measured, controlled and 
zeroed. We have presented a simple two-particle model 
that models the density shift data qualitatively. 
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APPENDIX A: SIDEBAND LINESHAPE 
DERIVATION 

The derivation in this appendix gene ralizes arguments 
presented by Boyd fld\, Le Targat [111, and Ludlow ^ 
to include longitudinal anharmonicity and population in 
higher bands. 

An atom's radial motional state (n^, Uy) shifts the side- 
band transition frequency via the r^z^ coupling term in 
the site potential. For a fixed longitudinal quantum num- 
ber Uz and detuning 5, only atoms with specific longitu- 
dinal quantum numbers (rix,ny) will be resonant with 
the probe laser. This resonance condition is given by the 
longitudinal energy gap 7 (Eqn. [5]). The sideband am- 
plitude at this detuning is determined by the Boltzmann 
factors for the corresponding quantum numbers. The 
single-TT-z sideband lineshape is a superposition of indi- 
vidual atomic transitions between states \i) = \nx,ny) 
and I/) = \nx,ny,nz + 1) , weighted by the Boltzmann 
distribution. The individual atomic transition is power 
broadened with base width F, given by the carrier Rabi 
frequency, since the natural transition linewidth is negli- 
gible. We also neglect higher-order contributions to the 
sideband where the change in longitudinal quantum num- 
ber is accompanied by a change in the radial quantum 
numbers since their amplitude is suppressed by an addi- 
tional factor of 7]^. We find 



'blue 



1 



linz)Y 



(Al) 



where Zr = exp[— ft,i^r/(^B7r)] is the Boltzmann factor 
associated with transverse state (riy) and transverse 
temperature Tr- By introducing a radial quantum num- 
ber Ur = Ux + riy and corresponding two-dimensional 
degeneracy factor (n^ -1-1), we can reduce the summation 
to one dimension: 



'blue 



(5)(X 



(rir + l)z^ 



rir + l 



^1 + ^[S- 7(n^) + :/rec^K + 1)] 



2 ' 



(A2) 

with base longitudinal gap 7(n^) = — Vrednz + 1). 

For detunings S/j{nz) <C 1 on the shallow sideband 
slope facing the carrier, the sum's main term will be as- 
sociated with the quantum number that minimizes the 
denominator, leading to a relation between detuning and 
quantum number: 



rir + 1 



-[S-liriz)]. 



(A3) 
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We approximate the sum by this main term and find 
6 



1 



e-"[i-*/^("^)le[7K)-(5], (A4) 



with a = slliil ^^j^ , \Ye ensure the appUcabiUty of the 
approximation by cutting off the hneshape at S = 7(^2) 
with the Heaviside function Q. For the relevant param- 
eter ranges, the hneshape can be approximately area- 
normalized to three significant figures with a prefactor 

The longitudinal trap anharmonicity places each of 
these single-band component lineshapes at slightly dif- 
ferent detunings. We find the final blue sideband line- 
shape as a Boltzmann-weighted superposition of area- 
normalized components: 



Here, the quartic approximation to the longitudinal po- 
tential reproduces the lowest energies and gaps to within 
a few percent. A better approximation to the energy En 
of the n-th longitudinal state is given by En / {hv-ccc) — 
[a„(g)-|-5„+i((7)]/2-f 2(7, where a„ and 6„+i are the char- 
acteristic values of the Mathieu equation boundingits n- 
th stability region at parameter value q (see Ref. 28] , for 
example). Here q has to be optimized numerically - start- 
ing from the harmonic approximation qo — Uq/ (Ahiy^cc) ~ 
such that the lowest gap {Ei — EQ)/h ~ j/i.cc[ai((7) — ao((7)] 
is equal to the measured longitudinal trap frequency. 



E 

(A5) 



Enjh = 1^2(^2 + ^) - -^("2 + "2 + !)■ 
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